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PACS 42 . 25 . Dd - Wave propagation in random media 

PACS 72. 15. Rn - Localization effects (Anderson or weak localization) 

PACS 72. 10. -d - Theory of electronic transport; scattering mechanisms 

PACS 72. 10. Bg - General formulation of transport theory 

Abstract - An amazingly simple model of correlated disorder is a one-dimensional chain of n 
potential steps with a fixed width l c and random heights. A theoretical analysis of the average 
transmission coefficient and Landauer resistance as functions of n and 5 = kl c predicts two distinct 
regimes of behavior, one marked by extreme sensitivity and the other associated with exponential 
behavior of the resistance. The sensitivity arises in n and <5 for 8 ~ tt, where the system is 
nearly transparent. Numerical simulations match the predictions well, and they suggest a strong 
motivation for experimental study. 
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Introduction. — An enormous amount is known 
about wave transport in disordered systems described by 
a potential with uncorrclatcd disorder (see, e.g., Refs. [HE] 
and references therein). The extension to studies of corre- 
lated disorder originated a couple of decades ago, when 
a number of surprising properties were found in one- 
dimensional (ID) systems. Correlated disorder studies 
have since multiplied considerably, and cover a broad 
range of topics. For example, in the random-dimcr model 
-a paradigm for the description of short-range correla- 
tions in ID disordered systems- fully transparent (de- 
localized) states have been discovered [3H6]. For long- 
range correlated disorder, a mobility edge in ID random 
potentials was found [7], and the authors construct po- 
tentials with mobility edges at specific energies; one can 
thus build systems having the desired transport properties 
by controlling the correlations. Perhaps more surprising, 
an Anderson-like metal-insulator transition [8] and viola- 
tions of single-parameter scaling due to short-range cor- 
relations [§] have been found. Rcf. [!T] developed Fokker- 
Planck equations similar to the ones found in the present 
paper. Dclocalization in the continuous random-dimcr 
model |10j . and in continuous disordered systems consist- 
ing of ^-potentials and barrier-well sequences [Hill] have 
also been discussed. 

Here, wave transport in ID disordered systems consist- 



ing of spatially extended scatterers -barriers and wells 
with a finite width- is considered, in contrast to previ- 
ous studies by our group [12] in which the spatial exten- 
sion of the scatterers and their separation were negligible 
and played no role in the analysis. The system contains 
n steps, assumed to be low with respect to the energy E 



< l c > 



< lc 



(r-l)/ c rl c 



Fig. 1: Schematic representation of an array of n steps of ran- 
dom height V r (r = 1, ■ ■ ■ n) possessing fixed spatial width l c . The 
incident energy E is taken larger than all the \V r \ 's. 

(see Fig. [T]) and characterized by: i) a fixed width l c which 
may fit an arbitrary number of wavelengths 5/2tt, where 
S = kl c , with k the wave number, and ii) random heights 
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V r (r = 1, ■■■n). The n heights V r are statistically in- 
dependent of one another, and identically and uniformly 
distributed, with zero average. In this model the corre- 
lations could scarcely be simpler, in the sense that for 
each step the potential is perfectly correlated within l c , 
but perfectly uncorrelated otherwise. One motivation for 
studying this model is that it could be interpreted (see 
Fig. [J) as simulating a potential described by a random 
process with a correlation length l c . Another motivation 
is that it exhibits intriguing transport properties which we 
now describe. 

Results and discussion. In the study of the model de- 
scribed above we find two distinct regimes of behavior. In 
what we shall call regime A, the gross-structure behavior 
of the average transmission coefficient (T) as a function of 
S and for a fixed number of scatterers n shows "bumps" 
near 5 = n, 2tt, • • • , and, correspondingly, the gross struc- 
ture of the average Landauer resistance [13] (R/T) (R= 
reflection coefficient; T= transmission coefficient) shows 
"valleys". This is illustrated in Fig. [2] for (T) and in Fig. 
[3] for (R/T), which include the region S ~ tt, this being 
the only multiple of tt that we shall consider in what fol- 
lows. For the case of weak scatterers, the system is almost 
transparent in regime A, and regime B (shown in Figs. [2] 
and [3]) is more localized. This gross-structure behavior 



ft, 




Fig. 2: Theory and numerical simulations for (T) vs 8, for a 

chain of n = 5000 scatterers, 10 5 realizations and for yo = 0.09 (the 
parameter related to the average step strength: see below Eq. (0). 
The main figure shows the gross-structure behavior as well as regimes 
A and B described in the text. The insets show the fine structure 
for 5 sa 7r and its enhanced behavior exactly at 5 = tt, leading to a 
dip at this point. This is consistent with the peak observed in Fig. 
[3] for the average resistance. The agreement between simulation and 
theory is excellent. The statistical error bar for 8 = n is ~ 10 — 5 . 

is not entirely surprising. The transmission coefficient T 
for a single barrier with fixed width and strength becomes 
completely transparent (T = 1) at the resonance values 
kl c = mr, n = 1, 2, • • • , where k is the wave number in 
the region of the barrier (5 > tt for low barriers). For a 
well, T = 1 at 8 < tt. For fixed step width and random 
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Fig. 3: Theory and simulations for the average Landauer resis- 
tance (R/T) vs (5 for the same system as in Fig. [2] The main figure 
shows the gross-structure behavior using a semilog scale, as well as 
regimes A and B. The inset shows the fine structure for 8 ss n with 
a linear scale. At 8 = tt is a peak, consistent with the dip observed 
in Fig. [2] for (T). The agreement between simulation and theory is 
excellent. The statistical error bar is the same as in Fig. [2] 



strength with zero average, (T) reaches a maximum value 
smaller than unity at 6 = tt. As the number of scatter- 
ers n increases, the "giant resonance" seen in the gross 
structure of (T) as a function of 5 is still similar to the 
above description for one random scatterer. However, for 
large n the system response becomes richer for values of 
5 very close to tt (this is a small subregion of regime A, 
to be denoted by S m tt), where it shows a remarkable 
fine structure. This is illustrated in Figs. [2] and 03 This 
fine structure is enhanced exactly at S = tt, where (T) 
shows a dip (Fig. [2|) and (R/T) a maximum (Fig. [3|); 
i.e., the trend of the system to dclocalize as it approaches 
S = tt from both sides reverses in an extremely narrow 
window around 5 = tt, where the system is less delocal- 
ized. The behavior of (R/T) can also be analyzed as a 
function of n for fixed 5. For regime B, 1 < 6 < 2.9 and 
S > 3.4, the familiar exponential increase of (R/T) with n 
is found;the mean-free-path (mfp) increases as regime A is 
approached, although this trend is reversed in the neigh- 
borhood of S = tt, consistent with the above description 
in relation with Figs. [2] and [3] The behavior of (R/T) as 
a function of n is shown for four values of S w tt in Fig. 
2] In this regime a small change in S changes drastically 
the behavior of (R/T) as a function of n: from a mono- 
tonic increase (<5 = tt) to an oscillating one (5 ~ tt), the 
wavelength decreasing away from 5 = tt. 

It is important to mention that gross and fine structures 
have been seen in the response of a variety of physical 
systems, and have been associated with various physical 
mechanisms. One of the first examples came from nuclear 
physics, where the structures were related to the many- 
body nature of the problem and gave rise to different time 
scales in the response [14]. In the present case, although 
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Fig. 4: Simulations for (R/T) vs n, using an ensemble of 10 

realizations, with yo = 0.09, and for four values of S ps it. (a) 8 = n. 
(b-d) 5 = 3.1405, 3.1400, 3.1380 (the results are symmetric around 
8 = it in the vicinity of this value). The statistical error bar is in 
all cases smaller than 10~ 6 and is not indicated. Also shown is the 
analytical solution of the differential Eq. I|13|l for the four values of 
8. 



a very precise analytical description of the numerically 
observed i) gross structure, ii) fine structure and iii) en- 
hanced behavior of the latter at S = -k is given ahead, 
an intuitive, physical explanation of ii) and iii) remains 
to be found. Indeed, Figs. [20 show the comparison of 
simulations with the theory. For the average resistance, 
the theoretical results are based on Eqs. ©, (fTTj) and 
the analytical solution of (|13|) ; for the average transmis- 
sion coefficient, on Eq. (|21[) and the explanation around 
Eq. (|22"j) . in regimes B and A, respectively. The agree- 
ment between the two is generally found to be excellent. 
This understanding and the extremely sensitive response 
suggest that experimental realizations of the model could 
lead to important devices. 

The theoretical model. — The theoretical treat- 
ment follows a familiar path. For the r-th scatterer (see 
Fig- P): define the dimensionless parameter y r = U r l 2 as 
a measure of its strength, where U r = 2mV r /h 2 . Since 
y r = (U r /k 2 )8 2 (recall that <5 = kl c ), for a repulsive bar- 
rier, U r > 0, 5 > y/j/7 means impinging above the top 
of the barrier. The transfer matrix for the r-th scatterer 
shown in Fig. [T] has the structure [2] 



a r fi r 

p; a* 



M r = 
due to time-reversal invariance. One finds 
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a r 
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2S 2 - y r 
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MV$ 2 ~ Vr). (2b) 



Landauer's resistance for this scatterer is R r /T r = \/3 r \ 2 - 

Landauer resistance of the chain. Consider next a 
chain of n scatterers to which one more, to be called a 
building block (BB), is added. Combining the correspond- 
ing transfer matrices, the average resistance for the chain, 
(l/3*- n ' ) | 2 ), obeys a recursion relation that couples with the 
quantity (qW/^W) as 



2<|/W| 

v(™)\ 



■)A(n) 



A(n + 1) - A(n) 
+2 [(a„ +1 /3* +1 )(/3Ha^) + c.c.j , (3a) 

( (n+i)£(n+i)} _ (a W^)) = (a n+1 p n+1 )A(n) 
+ ((al +1 ) 1) («(»>/*<»>) + (/3 ? 2 l+1 )(a(«)/3("))* . (3b) 



We have defined 



A(n) = l + 2(\/3^\ 2 } 



(i) 



where (• • • ) indicates an ensemble average and c.c. denotes 
the complex conjugate. The nature of the ensemble fol- 
lows from the definition of the statistical model: each y r 
is uniformly distributed over the interval (— 2/0)2/o) (then 
(Vr) = 0, (y 2 ) = 2/ 2 /3) and different y r 's are statistically 
independent. The recursion relations, Eq. ([3]), are exact 
and thus take into account all multiple scattering processes 
occurring in the chain. 

A comment on the scaling property of the model param- 
eters is in order. Although the transfer matrix of a single 
scatterer r depends, in principle, on the 3 parameters E, 
U r and l c , Eq. ^j shows that these 3 parameters occur in 
the combinations S and y r - Thus, for a specific realization 
of disorder, the resistance (R/T)^ of an n-step chain is a 
function of the n + 2 parameters 8,n,yi,---y n . An average 
over realizations, i.e., over j/i, • • • y n , is performed with the 
distribution function just described, which depends on the 
parameter j/o only. Thus the resulting average resistance 
((R/T)^) depends only upon the combination of 3 pa- 
rameters: S, n and j/o = Uq/1 2 , instead of 4. In Figs. [5] 
and [31 j/o and n are kept fixed and S is varied. This could 
be realized, e. g., by using a fixed l c throughout the chain 
of fixed length n and varying the energy. In each of the 
panels in Fig. [4j j/o and S are kept fixed and n is varied. 
This could be realized, e. g., by fixing l c and the energy, 
and varying the length of the chain n. 

Equation (J3]) admits two distinct approximations, which 
aid analytical treatment and correspond to regimes B and 
A, respectively. Define 



A'i 



L /3* +1 ) (/3("V">) 



(1/3, 



[l + 2(|/?(«)p)] 



(5) 



which is the ratio of the "coupling" containing (fi^a^) 
to the "direct part" containing (\P^ n ' | 2 )) on the right-hand 
side of Eq. (|3aj) . Numerically, it appears that the coupling 
can be neglected in regime B: for n < 5000 (larger values of 
n have not been checked for the analysis of K\) and for the 
value yo = 0.09 (corresponding to weak scatterers) to be 
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used throughout this paper (except in Fig. [5]), K\ < 1%. 
Dropping the coupling terms on the RHS of Eq. (f3"aj) gives 



A(n + l)-A(n) = 2(\f3 1 \ 2 )A(n), 



(6) 



since (|/3„ + i| 2 ) = (|/?i| 2 ). The solution of the recursion 
relation ((6j) with the initial condition A(0) = 1 is 

A(n) = (1 + 2(\f3 1 \ 2 )) n = e 2n i K^l*!')) = e 2 «W*. 

(7) 
This is the well known exponential increase found by Lan- 
dauer [13], with 



| = iln(l + 2(|/3 1 | 2 )) 



(8) 



For (|/3i | 2 ) = (Ri/Ti) <C 1, l/£ is approximately the av- 
erage reflection coefficient per unit length, which is iden- 
tified with the inverse mean free path (mfp) [12| and, in 
the present ID problem, is of the order of the inverse lo- 
calization length. 
Landauer's average resistance for the chain is thus 



(|/?(")| 2 ) = I(e 2 "^-l). 
Using Eq. ([2b"]). the quantity (|/?i| 2 ) is 

(l/3i| 2 H 



III 



AS 2 (6 2 - yi ) 



sin 2 y^ 3 



2/1 



(9) 
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where the ensemble average has to be taken with the distri- 
bution of yi described above. It can be evaluated exactly, 
but for small values of yo/5 2 we find 



(Ifl 



I, 



_!£ + O f|0) !£ = i«l 



f 



J 2 



Z c _ y 2 sin 2 (5 



£ 12 (5 4 



(11) 



where the approximate expression £ for the mfp shows ex- 
plicitly its dependence on the strength yo and the spatial 
extension l c of the steps. For the value yo = 0.09 used 
here, the approximate expression (fTTj) is an excellent ap- 
proximation to the exact result. 

The result of Eq. ^ , with the approximate mfp given in 
Eq. (fTTj) . was used to plot the theoretical average Landauer 
resistance in Fig. [3] for 1 < 5 < 2.9 and 5 > 3.4 (regime 
B) and is seen to give an excellent description of the data. 
Indeed, the results of the numerical simulations literally 
lie on top of the theoretical ones. Therefore regime B is 
marked by the lack of importance of the coupling terms 
(our first approximation). 

It is not possible to treat the behavior in regime A 
(2.9 < S < 3.4) without accounting for the coupling terms 
in Eq. ([3]). In particular, well inside regime A (5 « n) one 
finds considerable fine structure (see Fig. |4} which cannot 
be explained without the coupling. Figure |4ji shows the 
average Landauer resistance precisely for S = ir as a func- 
tion of the number of scatterers n, obtained from a sim- 
ulation. The approximate inverse mfp gives l c /£ = for 



5 = ir, and the theoretical result of Eq. © with £ replaced 
by £ would give (I/?*-™-*! 2 ) = 0, in gross contradiction with 
the increasing behavior with n obtained numerically. The 
inverse mfp of Eq. ([5]) using the exact average, Eq. (p~0|) . 
is non-zero for 5 = n (it is found to be of O(10~ 9 )), but 
not large enough to explain the numerical data. Indeed, 
for S in the edges of regime A, the coupling terms in Eq. 
(|3aj) become important and, for S sa n, dominate and the 
uncoupled difference equation §§§ becomes a bad approx- 
imation (also, it does not reproduce the result shown in 
the inset in Fig. [3]). 

This motivates the second approximation, which is a 
continuum approximation that converts the recursion rela- 
tions, Eq. ([3]), into differential equations. However, to do 
this the ensemble averages {\fi^\ 2 ) and (/3(™W n ') must 
behave "smoothly" as functions of n. For this purpose 
consider the recursion relations similar to those of Eq. ([3]), 
except with a BB containing m scatterers. It is their be- 
havior as a function of m which is of interest. Evaluating 
analytically the various BB quantities appearing in these 
equations in the same approximation as Eq. (fTTj) . i.e. ne- 
glecting terms of order (myo/S 2 ) 4 , leads some of them to 
have the form g(5)(e 2lmS — 1). For example, we find 



/ a* \ mic f (x \ 2in6 , nf m y°Y 

(a n+hn+m /3 n+hn+m ) = -y/i(5,m)e +°{-p-) - 

(12a) 
sin 5 - Se lS e 2imS - 1 



fi(5,m) = iS 



sin 3 S 2imS 



(12b) 



This result is, in general, not linear in m and thus pre- 
vents us from making a continuous approximation to con- 
vert the difference equations into differential equations. 
However, for S — it these BB expressions either van- 
ish or are proportional to m, and for S ~ n (regime A) 
they are all proportional to m. In these cases we have 
also verified numerically that (|/3( n+1 )| 2 ) w (|/3 (n) | 2 ) and 
( / g(™+i) Q ,(«+i)) ra (£(«) a («)) ) i.e., that these quantities be- 
have smoothly as functions of n, when n is changed by one 
unit. Thus, in this regime it is justified to make a contin- 
uous approximation, which gives the coupled differential 
equations 



3A{n) 

On 
db{n) 



A(n) + h(5)b(n) + f?(5)b*(n) (13a) 



On 



(e 



2M 



= 2 4 
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- l jf 4 (5y iS A(n) 
- l jh{S)e™ 



n) - ih(5)e M b*(n) 



(13b) 



neglecting terms 0(?/o/<5 2 ) 4 . Here, fi(6) is obtained from 
the above fi(5,m) setting m = 1, and fa(5), /3(d), fi{S) 
are similar functions of 6. The quantity A(n) is defined in 
Eq. 0]) and 

b{n)=e 2m5 {a {n) [3 (n) ). (14) 
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Equation f| 13[) is subject to the initial conditions 
A(0) = 1, 6(0) = 0. 



(15) 



The explicit expression for f^iS), f 3(8), f 4(8) and the exact 
analytical solution of Eqs. (fT3")) will be given elsewhere. 

That exact solution was used to plot the theoretical av- 
erage resistance as a function of 8 in Fig. [31 in regime 
A (2.9 < 5 < 3.4). Since the coupling has been taken 
fully into account, the agreement between theory and sim- 
ulations is excellent. Thus Eq. ([T3| gives a theoretical 
quantitative description of the fine structure, including the 
prominent peak at 8 = 7r shown in the figure. 

The analytical solution of Eq. (fT3|) was also used to plot 
(R/T) as a function of n in Fig. S^i-d. The agreement 
with the simulations is excellent exactly for 6 — 7r; the 
remarkable oscillatory behavior for 8 sa -k is a result of 
the coupling, and is generally reproduced quite well by 
the theory; the agreement deteriorates further away from 
8 = -k (see also right inset in Fig. [2]). There is evidence 
that the deterioration starts for larger n's as yo decreases. 
Agreement is also excellent for 5ft (a/3) ( n ) and 3(a/3)( n ) for 
1 < n < 5000 and for the same values of 8 shown in Fig. 
0] (these results will be presented elsewhere). 

Transmission coefficient of the chain. Finally, we an- 
alyze the transmission coefficient behavior for the chains 
studied above. It is convenient to use the polar represen- 
tation of Ref. [2] and introduce the notation A r = |/3 r | 2 for 
the r-th scatterer and A^") = |/3(")| 2 for the chain consist- 
ing of n-scatterers. In this notation, Eq. © is written as 
[seeEq. JTTJ] 



A(A< 



n)\ 



l + 2(A< n >>, 



where 



L 



AS = (Ai) = | 

Its continuous approximation is 
d(X) s 



ds 



1 + 2(A), 



(16a) 

(16b) 
(16c) 

(17) 



The evolution of (A), coincides with that found from the 
evolution equation for the A-probability density, iuj(A), 
known as Melnikov's equation [2] 



dw s (X) 
ds 



d_ 



A(l + A) 



dw s (\) 



OX 



(18) 



A similar result holds true for the second moment of A; 
a quantity K2 can be defined, analogous to K\ defined in 
Eq. ([5]), which is also found numerically to be small (< 
1.6%) in regime B. It is not possible to verify numerically 
the values of the analogous quantities K p (defined for the 



p-th moment) for all moments; but in fact it is not possible 
even for an individual moment if the order p is too large, 
because numerical control is lost in the calculation due to 
the rapid exponential increase of these moments with s. 
We thus propose the approximate validity of Melnikov's 
equation for regime B; this assumption allows finding the 
statistical properties of T which, in terms of A, can be 
written as 

t = TTa- < 19 > 

From Melnikov's Eq. (fl"8| . the expression for the p-th 
moment of the transmission coefficient can be reduced to 
quadratures, with the result |15j 



(T p ) = 



2c- s / 4 



r(p) Jo 

which, for the first moment, gives 



T(p - i + it) 



t tanh(7rf)d£, 
(20) 



(T) = 2e"^ 4 / 
Jo 



'7rt[tanh(7rf)/cosh(7rf)]dt. (21) 



This result was compared with numerical experiments in 
Fig. [H for 2.5 < 8 < 2.9 and 3.4 < 8 < 4 (regime B). The 
agreement is excellent. 

In regime A, the theoretical analysis uses 



(7>1-(A), 



(22) 



since (A) <C 1 (Fig. [3]), and (A) is taken from the analytical 
solution of Eq. (fl~3)) . This approximation for (T) was 
compared with the numerical simulations in Fig. [2] The 
agreement is excellent. The pronounced dip observed for 
(T) at S = it is consistent with the peak observed in Fig. 
[3] for the average resistance. 

To complement the above analysis, we may consider 
chains in which the potential strengths V r are not dis- 
tributed symmetrically around zero, and thus have a non- 
zero average. When all the V r > 0, one expects the 
"bump" in (T) of Fig. [5] to move to S > w, and, when 
all the V r < 0, to 8 < ir. This behavior is indeed found in 
numerical simulations, which also show a similar shift in 
the position of the dip, along with a loss of the symmetry 
around this shifted position. 

Finally, a few comments are in order with regards to 
the asymptotic behavior of our theory for large n's and 
the validity of the theoretical analysis. 

a) If the calculations of Figs. [2] and [3] are repeated for 
various n's: i) The "bump" in (T) and the "valley" in 
(R/T) become narrower as n increases from n = 1000 to 
n = 5000; it has been found, albeit with a smaller num- 
ber of realizations (1000), that this tendency persists to 
n = 50000. ii) The maximum in the bump in (T) remains 
undisturbed, iii) The peak in the average resistance for 
8 = 7r becomes ever higher as n is increased, in full ac- 
cordance with the behavior shown in Fig. [4ji, while the 
dip in the average transmission for 8 — tt becomes ever 
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deeper as the chain is made longer. Effects i), ii), hi) are 
reproduced well by the theory. 

b) As for the fine structure oscillations seen in Fig. [3] for 
{R/T) as function of 8 for a fixed number n = 5000 of scat- 
terers, the calculations have been extended to n = 20000. 
As n is increased, various features are observed (zoom- 
ing in on the results, to make them visible) in the region 
away from 8 = it where one can define a wavelength and 
an amplitude of the oscillations (for n = 5000 this region 
would be outside the extremely narrow window shown in 
the inset in Fig. 3): i) the wavelength decreases; ii) the 
amplitude decreases; hi) the extension of the region away 
from 8 = it where the oscillations are present before damp- 
ing out, decreases. An example of these three features is 
shown in Fig. [5l where an excellent agreement between 
theory and simulations is also seen. 



tion: a fine structure behavior which is enhanced exactly 
at 8 = 7r, where the system becomes less delocalized. In 
this region, a small change in 8 modifies drastically the 
behavior of the average resistance as a function of n. All 
of these phenomena are seen in simulations and are de- 
scribed very well by the theoretical analysis. We want to 
stress that our theory is a fully analytical theory, which, in 
regime B, is given by decoupling Eqs. ((3j), and in regime 
A by the coupled Eqs. (flU]) which, in turn, represent a 
continuum approximation to Eqs. ([3]). The theory has 
full predictive power, with no adjustable parameters. 

The success of our theoretical analysis and the extreme 
sensitivity we have described suggest the importance of 
the system experimental realization. 
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Fig. 5: Theory and simulations for the average Landauer resis- 
tance (R/T) vs <5 for three chains, with n = 5000, 10000 and 20000 
scattcrers, respectively, and 10 5 realizations. The decrease, with in- 
creasing n, of the wavelength and amplitude of the oscillations and 
the region where they are conspicuous, is clearly seen. The value 
3/0 = 0.01 was chosen to show these features more clearly. The 
agreement between theory and simulations is excellent. The statis- 
tical error bar is < 10 -9 and is not indicated. 

c) The fact, apparent from Fig. [4j that, as function of 
n (for fixed 8), the fine structure oscillations damp out as 
n increases, has been seen to persist to n — 100000 and is 
reproduced well by the theory. 

d) It is to be remarked that our theoretical analysis is 
restricted to small values of yo, as was indicated right after 
Eq. 1|13|). Therefore, although one expects, for a fixed yo, 
theory and numerical simulations to deviate starting from 
some Um-i there is evidence that Um can be made larger 
by decreasing yo- This fact was already remarked at the 
end of the previous subsection. 

Conclusions. — To summarize, we have discussed the 
problem of wave transport in ID disordered systems con- 
sisting of barriers and wells with a finite, constant width 
l c , and random strength. For weak scatterers, the system 
is almost transparent for 8 = k£ c ~ tt, and less delocalized 
farther away. For 8 ~ n, one observes a remarkable situa- 
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